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CHAPTER I 
INTRODUCTION 

In the recent past, several investigators have 
taken under consideration the phenomena observed when 
Charged particles interact with crystal lattices. 
Anomalous energy loss phenomena were first noticed in 
experiments on radiation damage in metals (1). Later 
reports cover several phenomena concerning the orient- 
ational dependence of reactions involving crystal lat- 
tices, including: 

1) Non-uniform and orientationally dependent 
penetration of thick crystals (2,3); 

2) Sharp orientationally dependent minima in 
large angle Rutherford scattering through thick crystals (4); 

3) Dramatic departures from predicted nuclear re- 
action yields in thick crystals (5); 

4) Non-smooth orientationally dependent angular 
distributions through thin crystals (6); 

5) Anomalously large orientationally dependent 
energy losses in thin crystals (7); and 

6) Anomalously low orientationally dependent 
energy losses of charged particles in thin crystals (8). 

While all these phenomena can be roughly categorized 
as channelling effects (9), the use of the term "Channelling" 





2 
has generally been restricted to the phenomenon of anon- 
alously low energy losses of charged particles in thin 
crystals. 

Charged particle channelling effects have been ob- 
served in silicon, germanium, and gold using protons 
(6,10,11), in germanium using deuterons (7), in gold 
using bromine and iodine ions (8), in aluminum using 
krypton ions and protons (2,3,11), and in copper using 
protons and krypton ions (11,12). In addition the slow- 
ing down of copper atoms in various lattice structures 
has been investigated by high-speed digital computer 
techniques (13). 

The investigations forming the basis for this thesis 
were made using thin monocrystalline films of mica, 
Silicon, and aluminum as targets for protons. The 
aluminum crystal was unsatisfactory, due to excessive 
thickness. Satisfactory results were obtained with the 
Silicon crystal, although due to equipment limitations 
the effect observed was not as dramatic as has been 
elsewhere reported (9). The primary result of the in- 
vestigations was the original observation of channelling 
in a multi-elemental crystal, mica. It is felt that 
this observation is of significance in the investigation 
of channelling. 








CHAPTER II 
THEORY OF CHANNELLING 

On entering a crystal along a direction of lattice 
Symmetry (row or plane), energetic positively charged 
particles experience a coulomb repulsion from the atoms 
at the lattice sites. For critical impact parameters 
relative to the center axes of the arrays, the part- 
icles are scattered through such angles that they be- 
come collimated along those directions in the crystal 
where 1) repulsion forces from arrays on either side of 
the particle path are in approximate equilibrium, and 
2) the electron density is relatively low. The traject- 
ories of all other particles become random within a few 
collisions. 

The theory of channelling is concerned basically 
with three predictions: 

1) prediction of channelling orientations; 

2) prediction of the limiting incidence angle 
for observable channelling; and 

3) prediction of the magnitude of the energy 


anomaly. 





“Due to the embryonic nature of channelling theory, 
no Single authoritative source work exists. No attempt 
is made, therefore, to fully reference this chapter. 
Primary sources for this chapter are references 14 and 15, 





Prediction of channelling orientations can be 
easily made through a basic knowledge of crystallography. 
Figure 1 shows the relationships between various planes 
and axeS in a typical hexagonal close-packed crystalline 
lattice. Both silicon and mica exhibit this type of 
structure. In mica the (100) axis is perpendicular to 
the face of the crystal, but no such easy identification 
is available for silicon. The assumed orientation for 
Silicon will be discussed in Chapter V. 

The prediction of the limiting incidence angle for 
stable channelling proceeds as follows: A charged part- 
icle incident at angle Y to a symmetry array in a crystal 
has components of momentum both along the symmetry array 
(longitudinal) and normal to it (transverse). During the 
interaction, for small ¥Y, the longitudinal momentum may 
be considered constant. The transverse momentum is mv‘, 
if Y is small, and the associated energy is BOY", where m 
and v are the mass and velocity-of-the incident particle 
and E, is its total incident energy. 

The transverse motion of a particle with this energy 
must be thought of as equivalent to a particle in an 
average potential field due to many atoms, for sufficiently 
large distances from the array. 

If Pain represents the distance of closest approach 
of the particle to a row, then it satisfied the PE pare 
Ship | 


a ae ey 
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2 Dieapee Oi eis 
1-(u LO ae yD eee) /E 4 ) = 0 (1) 


where V(P an) is defined as the average row potential at 
a distance Prin? and u is the impact parameter of the 
interaction in a hypothetical plane perpendicular to the 
row on which is projected the trajectory. 
If the partic trajectory is in the plane of the 

row of atoms, u = O and 

ee Vipe. > = 

O min 
Thus the limiting incidence angle for stable channelling 
is defined, if the closest distance of approach is known. 
Although (2) is only correct in the extreme situation 
Where u = O, it is exact for an atomic plane where Y 
is the angle of incidence of the particle to the plane. 
But defining act and determining ¥Y, depend on the 
proper choice of a potential form. Table I shows aver- 
age row and plane potentials for different interaction 
potentials V(r) between a charged particle and a single 
atom. 
The determination of the magnitude of the channelling 

effect involves the definition of an anisotropy coefficient, 
KF, as the fraction of the transmitted beam with anomalously 


low energy loss. An empirical formula for F is given by 


= 22 ep 
r= L-exp | -[(25/(m 0, )*) (exp(-¥ /0, ) 1] (3) 
Where ¥ is the angle of incidence to a planar channel, $ 


is the half-angle of acceptance of the channel, and 0° 





fi 00] 





Figure 1 





7 
is the mean-square angle of multiple scattering through 
thickness t (14). For a given value of t the fraction F 
varies with y approximately as exp(-¥°/8,*). Thus a slow 
decrease of F with increasing ¥ is to be expected. 


Table II compares values of F found experimentally 





here with values computed using equation (3). As can be 
seen, the comparisons indicate that equation (3) provides 
only a rough approximation of F. 

The problems in fitting this theory to experimental 
data have been generally in the choice of a potential 
form, allowances for thermal vibration of the crystal, 


and aberrations due to non-uniformity of the potential 


field. 
The screened coulomb potential of Bohr, 
V(r)=@ z,e°/r)exp(-r/a) (4) 
and of Born and Mayer, 
V(r)=Aexp (-r/a) (5) 


are most often used, but neither gives good results over 
the entire range of particles and targets used in experi- 
ments to date. The conclusion may be drawn that a single 
potential expression valid for all cases would involve 
Some combination of existing potential forms, with ap- 
propriate range form factors. 

The effects of thermal vibration depend on the dis- 
tribution of displacements about a lattice site. Assun- 


ing a Gaussian distribution with variance of ae normal 
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9 
to the particle trajectory, a Born-Mayer potential is 
transformed to 

V'(p)=V(p) exp (ur /2v) (6) 
where b is some constant, and a1° is a function of the 
Debye temperature of the target. It is seen that V'(p) 
increases slightly with increasing temperature, leading 
to a corresponding increase in ¥ (14). No temperature~ 
dependence studies have as yet appeared in the literature. 

The foregoing discussion indicates that the theory 
of channelling is quite inadequate. 





CHAPTER IIIl 
APPARATUS AND TECHNIQUES 


Equipment 


The proton beam used in this experiment was provided 
by a model KN-4000 Van de Graaff accelerator manufactured 
by the High Voltage Engineering Corporation. The accel- 
erator is capable of providing protons of from 0.5 to 
4,0 million electron volts energy at a maximum beam 
current of 400 microamperes. 

The accelerator is mounted vertically so as to 
produce the beam heading downward. Directly below the 
accelerator is a ninety degree momentum analyzing magnet 
Which deflects to the horizontal only those particles of 
@ momentum corresponding to the desired energy. 

This horizontal beam is focused down the beam tube 
by a quadrupole focusing magnet (See Figure 2). Two 
fixed tantalum collimators and an adjustable slit tant- 
alum collimator provide additional definition of the 
beam spot on the primary target in the scattering 
chamber. 

The scattering chamber is an aluminum cylinder, two 
feet in diameter and 10.5 inches high. The primary 


target holder is located vertically on the axis of the 
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Chamber. A target mounted on the target holder is thus 
directly in the path of the incoming beam. The vertical 
target may be rotated through 360° about a vertical axis 
without disturbing the system vacuum. 

In the sides of the scattering chamber are ten two- 
inch diameter ports, to permit entry and exit of the bean, 
to attach the vacuum pumps, to measure inside pressure, 
and to permit observation of the primary target through 
clear windows. Unused ports are covered with blank brass 
plates. In this experiment a Faraday cup was fitted on 
the straight-through port, the channelling apparatus to 
be described was placed on the 45° (from straight-through) 
port, and all other ports, with the exception of the beam 
inlet port, were unused. 

The Faraday cup collects the remnants of the beam 
passed through the primary target. The charge thus ac- 
cumulated is measured by a current integrator, and this 
measurement is used to calculate the total amount of 
beam passing through the target. The amount of beam 
lost in the target due to scattering and nuclear re- 
actions is negligible compared to the total beam because 
of the thinness of the primary target used. 

The channelling apparatus (See Figures 3 and 4) 
affixed to the 45°-port on the scattering chamber con- 
sists of four sections. In order outward from the 


scattering chamber they are: 
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1) <A ten-inch section of two-inch brass pipe; 

2) <A brass flange with a tantalum collimating 
disc mounted in its center; ~~~ © 

3) A crystal target holder assembly; and 

4) A detector and monitor assembly. 

The ten-inch section of two-inch brass pipe assists 
in the collimation of the secondary beam by increasing 
the base line for beam dispersion calculations from 12.5 
inches to 22.5 inches. 

The brass flange supports the first collimating 
disc. This disc is a thin sheet of tantalum with a 
small collimating hole in the center. Four discs can 
be interchanged to give a hole size as small as 0.013 
inch. In this experiment a 1/32-inch hole was found to 
be the smallest hole which would allow a sufficiently 
rapid counting rate. This size collimating hole 22.5 
inches from a 1/8-inch diameter beam spot on the primary 
target gives a half-angle of dispersion on the secondary 
target of Ont. 

The crystal holder assembly consists of a short 
length of two-inch brass pipe with a smaller brass tube 
inserted in its side normal to the axis. The crystal 
holder is mounted on a brass rod which moves inside 
this tube. A flange on the tube and a plate on the 
upper end of the rod have angular markings around their 


rims by which measurements of the angle 9 to a relative 
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15 
accuracy of 0.1° can be made. An O-ring on the rod allows 
adjustment of crystal position in 9 without disturbing the 
system vacuum. The end flanges on the pipe also have 
angular markings on their rims, whereby the angle yY can 
be measured to a relative accuracy of 0.1%. Here again, 
with O-ring seals at either end, » settings may be 
Changed without affecting the system vacuun. 

The detector and monitor assembly consists of an 
at the wide end, with two monitor detector ports mounted 
in the sides at the narrow end and a sliding plate over 
the open end. The two monitor detector ports are brass 
tubes placed at an angle of 25.2° with the beam axis, 
and so located that their axes converge at the target 
position when the crystal holder assembly is in place. 
Inside the brass tubes slide brass rods with provisions 
for mounting monitor detectors on their inner ends. The 
tubes are adjustable so that the detectors may be pos- 
itioned as close as one inch or as far as five inches 
from the target. O-ring seals on the brass rods allow 
repositioning of the monitor detectors without breaking 
the system vacuum. In the sliding plate at the wide 
end are mounted the detector collimating device and the 
Solid state detector. This plate allows the positioning 
of the detector in the X-Y directions normal to the beam 
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16 
path. The solid state detector is mounted in a recess 
in the center of the plate, with the collimating disc 
mounted one inch in front of it. This disc arrangement 
is identical in construction to the first collimating 
disc. Again it was found that a 1/32-inch hole was 
optimum. The plate is attached to the end of the as- 
sembly by four clamps, which are independently driven 
by screws. Two of the screws, as can be Seen in Figure 
5, are Calibrated, and turn along a calibrated rule. 

By these screws the position of the sliding plate, and 
thus the second collimator and the detector, can be 
determined to an accuracy of 0.001 inch in X and Y. 
The solid state surface barrier silicon detector (ORTEC 
Model SBEJO25) is mounted behind the collimator on the 
Sliding plate. 
Electronics 

The signals from the detector mounted on the slid- 
ing plate are fed into the low noise preamplifier 
section of an ORTEC Model 101 amplifie® (See Figure 5). 
After being amplified in the preamplifier the signals 
are passed through a biased amplifier, which further 
amplifies them. The resulting signals are analyzed in 
a Nuclear Data Corporation Model ND-110 128-channel 
pulse height analyzer. The resulting analysis is dis- 
played on a semi-logarithmic piot on a cathode ray tube, 


With number of counts as the ordinate and channel number 
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17 
as the abcissa. tach channel in the analyzer represents 
an energy Window in the spectrun. 

By varying the gain and the bias on the biased 
amplifier the portion of the energy spectrum to be 
analyzed and the width of the individual energy windows 


may be varied. 


Target Construction 

As previously mentioned, thin monocrystalline mica, 
aluminum, and silicon targets were used. 

Mica is readily obtainable in very thin, uniform 
Sheets. Four targets were used, one of which was a 
commercial mica window. The other three were prepared 
by picking apart a "book" of pressed phlogopite mica 
oe we ater With an ordinary straight pin and 
razor blade. The valuable feature of mica is that 
naturally occurring "books" of mica cleave along the 
base direction perpendicular to the (100) axis. 

The preparation of the aluminum target involved 
electropolishing a cast aluminum plate to the desired 
thinness. A functional drawing of the apparatus for 
electropolishing is shown in Figure 6. The aluminum 
plate is used as one electrode, and an aluminum foil 
as the other. The proper electrolytic solution for 
electropolishing aluminum consists of one part nitric 
acid and four parts methanol. Since this mixture is 


explosive at room temperatures, it must be cooled with 
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dry ice. A potential of approximately six volts D0 
between the electrodes produces the desired effect, 
that is, the removal of aluminum from the plate. Un- 
fortunately, successful electropolishing requires a 


great deal of practice, and time was not available to 


._  _ae = 


obtain the skills necessary to produce crystals of the 
thinness desired. 

Silicon targets are prepared by chemical etching. 
The solution required consists of 85 ml nitric acid, 
DeemienvydroLiueric acid, anda mis acctic acid. This 
Solution polishes at the rate of one mil per minute. 
Extreme care must be taken in the final stages of 
either electropolishing or chemical etching, because 
the thin crystals have very little strength, and crumble 
easily. | 
Zechniques of Operation 


A system vacuum better than 10 


2 


mm Hg, monitored 
With thermocouples and ion gauges, is required for the 
operation. Special venting ports in the first collin- 
ator flange and in the detector plate assembly facil- 
itate pumping out the apparatus. 
When the main vacuum isolating valves have been 

opened, and the particle beam is directed into the 
Scattering chamber, the beam spot size must be reduced 


to about 1/8 inch-diameter on the primary target. This 


a —_ 


—— a = == 








Aluminum 
Pm crystal 
Nitric Acid- 


Methanol Blectro- 
ae cooled by 
BB gt a 


=a 


Rheostat 





1 
! 
! 
! 
! 
' 
’ 
I 
! 
! 
’ 
' 
! 





Aluminum 
ee jails 


Figure 6 


20 





on 
is accomplished either by adjusting the movable slit 
COLL mowers or by chenging the field of the quadrupoxe 
focusing magnet. 

Once the beam spot is as desired with the channel- 
ling target out of the beam path, the slidince plate is 
moved in X and Y until a maximum count rate is regist- 
ered on a rate meter connected to the analyzer. Lach 
time the beam is turned on, this position must be re- 
located because it changes with the slightly varying 
characteristics of the beam and the primary beam spot. 
When the X-Y position of maximum count rate has been 
found, the beam axis has been geometrically defined. 

The taking of data involves first measuring an 
energy spectrum of particles reaching the detector 
without passing through the channelling target. By 
knowing the initial energy of such particles (from the 
settings on the momentum analyzing magnet), and their 
final energy, the energy loss in the first target can 
pe devernimoed., “Thiis determination facilitates the 
Coulci. Sumron Of the Crreetive thieékness of the first 
target. Once the effective thickness is known, the 
Sspeciite=eneray Wess ol any type of particles at any 
men demu cHemayewcon beeioundebywcalculatiton, or in 
most cases by consulting existing tables. 

It must be noted that effective thickness depends 


also on the target orientation with respect to the 
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incident beam. Care must be taken to reproduce the 
orientation of the first target to within a few degrees 
each time, to minimize fluctuations in data. 

After the relevant data without the second target 
have been obtained, the target is lowered into position, 
at an arbitrary angle of 9, where @ is the angle bet- 
ween the beam path and the normal to the crystal face. 
several runs must first be taken at various angles of 
9, to enable measurement of secondary target thickmess. 
farget thickness, in energy units, is simply 

we (E,-5,) cos @ (7) 
where Es is the energy of particles incident on the 


second target and &. is the energy of particles incid- 


d 
Emuwom unc Oeucceuor. FaGLOrTS to convert t from units 
of energy to units of ne/em® are readily available in 
tabular form (16). Correction of data for target thick- 
HeSseasean imporcant step in the analysis of data. 

fhe search for channelling in the target crystal 
involves taking gata in 0.1° or 0.2° steps about 2° 
either side of the predicted orientation. As previously 
neoted, a single channel analyzer and a scaler count the 
total number of particles striking the detector on each 
run, thus providing normalization. At the conclusion 
of each run, the spectrum is recorded on a numerical 
foie OU ana atco ON a curve tracer. The resulting 


curves are thus available for immediate comparison with 
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Ogee cuns.. 

initial runs recorded 1000 counts on the scaler, 
Tom a2 Stavisticai error of less than three percent. 
After the channelling orientation had been located, 
10,000-count runs were taken at the channelling angle 
and at an orientation just beyond the limiting angle 
for observable channelling, to give a statistical 


error of the order of one percent. 





CHAPTER IV 
TREATMENT OF DATA 


imeatmentwot datarinvolived three distinct processes; 

1) Selection of best raw data; 

2) Application of corrections to best data; and 

3) Selection of best display of results. 

The selection of best raw data was made by direct 
comparison or the crude curves drawn by the curve plotter 
at the end of each run. This method was found to be more 
facile than comparisons of numerical readouts. 

The data thus selected were corrected only for tar- 
get thickness. As relatively little information on 
peak coureddgine due to multipie scattering or energy 
Susepeling in Such cases is available, no attempt was 
made to correct for these effects. It is felt, however, 
that any contributions they might make would be small 
Cormarcasuoneje Channelling eifiects observed, and would 
remain nearly constant through the small ranges of 0 
over Which channelling was observed. 

Poa etornLNc Qaya to be useful the numerical 
Peep OULTecOunus per Channel must be converted To 
counts per energy increment. To accomplish this it was 


necessary to take a sequence of electronics-only runs, 
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25 
using as 2 pulse source a calibrated pulse generator 
built into the amplifier. A representative graph of 
data taken is shown in Figure 7. 

ties elt tnat 2 Stralghtiorwamd Camtesian plot 
of results would be most meaningful. A semi-logarithmic 
plot shows a more pronounced high-energy shoulder, but 
has tne effect of distorting the magnitude of the smail 


effects observed. 
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CHAPTER V 
RESULTS AND CONCLUSIONS 


Results 

Channelling was obsexzved in mica and in silicon. 
In mica the effect was observed on, or very near, the 
(100) axis normal to the face of the sample. In silicon 
the effect was observed on the (i111) plane at an angle 
oF about Beh to the (i100) axis. This determination was 
besed on the axis information obtained from the source 
Peco Crystal~za.tne department .of blectrical Engineering. 

In mica, the limiting angle for observable channel- 
ling was found to be 0.3° + 0.05°, and 9.2 percent of 
Hone particles emerging from the crystal at the channel- 
ling orientation were found to have anomalously high 
E@ewey. in silveon, the limiting angle for observable 
channelling was found to be 0.3° + 0.05°, and 6.8 percent 
of the emergent particles had anomalously high energy. 
Figures 8, 9, and 10 show these results. 

For both materials the anisotropy coefficient F 
was calculated (using equation (3)), and resulting values 
were compared with the experimentally obtained values. 
These values are shown in Table II. 

Cellevwationus Of tne limiting angle for observable 
channelling were made using equation (2). A comparison 
of these values with the actuel observations is shown in 
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Table If 


Anisotrovy coefficient (F) 


Material (predicted) F (observed) 





Table III 


Limiting incidence angle for observable channelling 


Material predicted observed 





* Predictions in mica are based on an average Z at each 
lattice site of 10.42, 
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oy 
Conuelusionus 

Pagires (, o,.ond 9 Lead to the conclusion that 
channelling has been observed. 

From Table II it can be seen that an empirical 
prediction of the anisotropy coefficient by equation 
(3) is only a rough approximation. 

fable ITI shows that the definition of the limiting 
angle for observable channelling by equation (2) is not 
exact. 

The final conclusion to be drawn is the one stated 
at the end of Chapter II: The theory of channelling is 
as yet too inadequate to safely predict tne magnitude 
or extent of the effect. Much experimental effort must 
be expended before sufficient data are available on 


which to base a workable theory. 
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